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Abstract. Muon spin rotation and relaxation (/xSR) experiments have yielded 
evidence that structural disorder is an important factor in many f-electron-based 
non-Fermi-liquid (NFL) systems. Disorder-driven mechanisms for NFL behaviour 
are suggested by the observed broad and strongly temperature-dependent fiSR 
(and NMR) linewidths in several NFL compounds and alloys. Local disorder- 
driven theories (Kondo disorder, Griffiths-McCoy singularity) are, however, 
not capable of describing the time-field scaling seen in muon spin relaxation 
experiments, which suggest cooperative and critical spin fluctuations rather than a 
distribution of local fluctuation rates. A strong empirical correlation is established 
between electronic disorder and slow spin fluctuations in NFL materials. 
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1. Introduction 

The observed breakdown of Landau's Fermi-liquid paradigm in a number of metallic 
systems has led to an explosion of effort to understand this non-Fermi liquid (NFL) 
behaviour [1-3]. The NFL breakdown is usually identified by deviations of low- 
temperature thermodynamic and transport properties [specific heat C(T), magnetic 
susceptibility x(T), electrical resistivity p{T)\ from the predictions of Landau Fermi- 
liquid theory [C(T) oc T, X (T) = const., Ap(T) = p(T) - p(0) oc T 2 ] [4]. Such 
deviations been observed in many f-electron heavy-fermion metals and alloys [2] , and 
often (but not always) take specific forms: the specific heat Sommerfeld coefficient 
7(T) = C(T)/T diverges logarithmically as T -> 0; x(T) varies as 1 - T 1 / 2 or - In T 
or T- rn (to w 1/3); and Ap(T) varies as T n ,n< 2. 

By definition no phase transition occurs at a nonzero temperature in an ideal 
NFL system. Nevertheless, in most NFL systems a magnetic phase is nearby in the 
appropriate phase diagram; the transition temperature T C (S) is reduced by varying 
an experimental parameter 6 (pressure, composition, magnetic field, . . . ), and NFL 
behaviour is observed in the neighborhood of the critical parameter value <5 C for which 
T c is just suppressed to zero. This is a quantum critical point (QCP) [5], which 
separates magnetically ordered and paramagnetic phases at zero temperature. Most 
theoretical effort in this area has concentrated on effects of the QCP. Magnetic 
resonance and other measurements have shown, however, that structural disorder is 
an important component of the NFL phenomenon in many systems. 

Nuclear magnetic resonance (NMR) and transverse-field muon spin rotation (TF- 
^iSR, i.e., magnetic field applied transverse to the muon spin) measurements in 
UCu4Pd and UCu3.sPdi.5 [6-8] demonstrated that the magnetic susceptibility in these 
NFL materials is strongly inhomogeneous. Following this discovery two broad classes 
of theories began to address the role of disorder in NFL behaviour. In the 'Kondo- 
disorder' approach [6,9-11] structural disorder gives rise to a broad distribution of 
Kondo temperatures Xk- NFL behaviour arises from low-Tk spins that are not in 
the Fermi-liquid state. In the 'Griffiths-McCoy singularity' picture [12, 13] spin-spin 
interactions freeze low-Tk f moments into clusters with a wide distribution of sizes; 
the larger clusters dominate the susceptibility and lead to divergent behaviour as 
the temperature is lowered. Both Kondo-disorder and Griffiths-McCoy singularity 
scenarios seem to be compatible with observed NMR and TF-pSR linewidths [14, 15]. 

Longitudinal- field muon spin relaxation (LF-/1SR) experiments in NFL 
UCus-xPds, x = 1.0 and 1.5 [16,17], and CePtSii-^Ge^, x = and 0.1 [18], 
indicate, however, that the f-ion spin dynamics in these alloys are better described 
by a picture in which critical slowing down occurs cooperatively throughout the 
sample, rather than at (rare) low- Ik spins or large clusters. The muon spin relaxation 
rates are nevertheless widely distributed, and the dynamic behaviour closely resembles 
that of spin glasses [19,20]. Thus an important question is whether disorder in the 
interactions dominates the dynamics, as in a spin glass above the glass temperature, 
or if, instead, the critical slowing down is still controlled by a QCP. Theoretical 
treatments exist [21,22] that combine elements of both these viewpoints, and describe 
a 'quantum spin glass' with a suppressed or vanishing glass temperature. 

In the meantime NFL systems have been discovered in which disorder does not 
appear to play an essential role. At ambient pressure these include CeCus.gAuo.i [23- 
25], Ce(Rui_ K Rh :c ) 2 Si2, x w 0.5 [26,27], CeNi 2 Ge 2 [28,29], and YbRh 2 Si 2 [30]. Thus 
one can compare spin dynamics in 'disordered' and 'ordered' materials. The result 
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of this comparison is evidence that the rapid muon spin relaxation in the disordered 
systems is due to slow dynamics associated with quantum spin-glass behaviour, rather 
than (homogeneous) quantum criticality. This evidence consists of 

• the observation of a wide inhomogcneous distribution of relaxation rates 
characteristic of disordered systems when the relaxation is strong, 

• the fact that in the ordered NFL compounds CeNi 2 Ge2 and YbRh 2 Si2, and even in 
the doped alloys CeCus.gAuo.i and Ce(Ruo.5Rho.5)2Si2, the low-frequency weight 
of the spin fluctuation spectrum is more than an order of magnitude smaller than 
in the disordered NFL systems, and 

• remarkably good correlation of the muon spin relaxation rate and its 
inhomogeneity with the residual resistivity of NFL systems. 

This paper reviews /iSR studies of NFL materials that shed light on the interplay 
between disorder and NFL behaviour. Candidate disorder-driven NFL mechanisms 
are discussed in section 2, which describes the way in which TF-^iSR (and NMR) 
linewidths probe the spatial inhomogeneity of the magnetic susceptibility that is 
the principal signature of disorder-driven NFL behaviour. As an example, TF-/iSR 
evidence for a disorder-driven NFL mechanism in the NFL compound UCu4Pd is 
described in section 2.2. Section 3 presents evidence that muon spin relaxation 
functions in at least some of these materials obey time-field scaling [19,20,31]: the 
dependence of the spatially-averaged muon spin relaxation function G(t, H) on time t 
and longitudinal field H obeys the scaling relation [19] 

G(t,H) =G(t/H 1 ) . (1.1) 

Time-field scaling is a signature of slow dynamics, i.e., a zero-frequency singularity 
in the fluctuation noise spectrum. It can in principle appear near any critical point, 
but is usually associated with spin-glass-like behaviour. In NFL materials such glassy 
spin dynamics may arise from the effect of disorder on quantum critical fluctuations. 
In section 4 muon spin relaxation results are compared with the residual electrical 
resistivity (a measure of disorder in the electronic system) and the low-temperature 
Sommerfeld coefficient (a measure of the characteristic electronic energy scale), and 
found to correlate well with the former but not the latter. These results establish 
strong links between quantum criticality, slow fluctuations, and electronic disorder 
in NFL materials [18,32]. The /xSR experiments and questions they raise raised are 
summarized in section 5. 

2. Disorder-driven NFL mechanisms 

As described briefly in the Introduction, NMR and TF-/1SR experiments [6-8] have 
revealed anomalously large and strongly temperature-dependent resonance linewidths 
in a number of NFL materials. It was realized [6] that this behaviour is characteristic 
of a spatially inhomogeneous distribution of local f-electron susceptibilities. This 
in turn suggests that a characteristic excitation energy A associated with the f 
electrons, which enters the zero-temperature susceptibility via the relation \ oc 1/A, 
is inhomogeneously distributed due to structural disorder. 

In the single-ion 'Kondo-disorder' picture [6, 9] this energy is the Kondo 
temperature T K , which is broadly distributed because of its extreme sensitivity to 
f-electron-conduction-electron hybridization. If the distribution function P(Tk) is 
broad enough so that P(Tk=0) ^ 0, as shown in figure 1, then at any nonzero 
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Figure 1. Distribution function P(Tk) in the Kondo-disordcr model. At a given 
temperature T impurities with Tk < T (shaded region) are not compensated 
and give rise to the non-Fermi-liquid behaviour of the system. This behaviour 
continues down to T = if P(Tk=0) is nonzero. After reference [33]. 



temperature T the f ions with Tk < T arc not Kondo compensated and are the 
source of NFL behaviour in the Kondo-disorder model. Bulk thermodynamic and 
transport properties are averages over -P(Tk) of the usual Kondo expressions for these 
quantities. It can be shown [9] that j(T) oc \{ T ) oc - InT, and Ap(T) oc T. 

A simple heuristic calculation may clarify how the inhomogeneous spread 5 X in 
X might be expected to behave in a Kondo-disorder model. Assume a spread STk in 
Kondo temperatures Tk. Assume further that the Curie- Weiss expression 

xCOocCT + Tk)" 1 , 



captures the essence of the local Kondo physics. Then 6 X IS given by 

(T + T K )- 2 5T K ot X 2 (T) 



(2.1) 

(2.2) 
(2.3) 

for small <5Tk- Thus b~x/x should be roughly proportional to X , with temperature an 
implicit parameter, reflecting the increasing effect of the distribution in Tk in equation 
(2.1) as the temperature is lowered. We do not expect this picture to be very accurate, 
as it assumes a small spread in Tk and a simplified Curie- Weiss susceptibility, but a 
strong increase in 5x1 X a t large x (l° w temperatures) is evidence for the basic effect. 

The Griffiths-McCoy singularity model [12, 13] also invokes a broad distribution 
of Tk, but assumes that uncompensated f-ion spins interact via RKKY coupling and 
form clusters. These clusters then tunnel between orientational configurations via 
Kondo spin-flip processes, and the distributed energy A characterizes these tunneling 
processes. The model exhibits Griffiths-McCoy singularities associated with a QCP 
below percolation threshold. The thermodynamic properties diverge as T — > in 
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a manner described by a nonuniversal exponent A < 1. Predictions for the specific 
heat C(T), the magnetic susceptibility x(T), and the fractional rms width Sx(T)/x(T) 
of the inhomogeneous susceptibility distribution are given by 

7 (T) cx X (T) cx T- 1+A ; (2.4) 

S X (T)/x(T) cx T- A /2 . (2.5) 

Thermodynamic properties have been shown to be in good agreement with this theory 
for a number of NFL alloy systems [34] . 

For both the Kondo-disorder and Griffiths-McCoy singularity models fits to the 
bulk susceptibility determine the distribution function P(x), from which Sx(T) / x(T) 
can be calculated. This can then be compared with NMR or TF-^iSR linewidth data 
with no further adjustable parameters [6-8, 12] as described below; good agreement in 
a given NFL material can be taken as evidence for applicability of a disorder-driven 
mechanism of this kind. This procedure often works equally well for both the Kondo- 
disorder and Griffiths-McCoy singularity scenarios [14, 15], so that other experiments 
are necessary to distinguish between them. 



2.1. NMR/TF-fiSR linewidth and susceptibility inhomogeneity 

We briefly review the relation between the static resonance linewidth and the spread 
in local susceptibility. More complete descriptions are given in references [8] and [33]. 

In a paramagnet in an applied magnetic field B the time-average or static 
local field (JBl) at the site of a spin probe (muon or nucleus) is shifted by an 
amount AB = (-Bl) — Bo due to the hyperfme coupling (direct or transferred) 
between the spin probe and the surrounding electronic spin polarization. The relative 
shift Ki = \AB\i/\B \ of the i th spin probe is given by 

Ki = ^2aijXj , (2.6) 

3 

where Xj is the susceptibility of the j th f ion and is the transferred hyperfme 
coupling constant between spin probe i and ion j. A spatially inhomogeneous 
distribution of the Xj gives rise to a corresponding distribution of the K i: and hence 
broadens the resonance line. 

From equation (2.6) the relation between the rms width Sx of the susceptibility 
distribution and the rms width SK of the shift distribution can be shown to be 

5K = a*S X , (2.7) 

where a* is an effective hyperfme coupling constant that depends on the spatial 
correlation of the disordered susceptibility and can be calculated from the [8] . Thus 
the spread in shifts is a direct measure of the spread in x, a quantity that is very hard 
to obtain from other experiments. Determination of a* in general requires knowledge 
of the spatial correlation function that describes disorder in the inhomogeneous 
susceptibility; this correlation function can be determined from spin-probe spectra 
only if data are available for more than one spin probe [7, 8] . The calculation simplifies 
in the extreme limits of long-range correlation (LRC) (correlation length 3> range 
of the hyperfme coupling) and short-range correlation (SRC) (correlation length < 
nearest-neighbor f-ion distance). It is found that 

a/2 



4 LRC 



Y.jOij and a| RC = (E^) . (2.8) 
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Thus one can find S\ from SK. Furthermore, the quantity 5K/(a*x) is an estimator 
of the fractional spread in susceptibility Sx/x- 

A convenient way to compare the data to these model results is to plot 5K/(a*x) 
[cf. equation (2.7)] and the theoretical $x(T)/x(T) (cf. section 2) versus x(T), with 
temperature an implicit parameter [6,8]. An example of such a comparison is given 
below for the NFL compound UCu4Pd. 

2.2. Evidence for a disorder- driven NFL mechanism in UCu^Pd [35] 

Heavy-fermion materials in the series UCu5_ x Pd K , 1 < x < 1.5, exhibit NFL 
behaviour [36,37]. These alloys crystallize in the fee AuBes structure (space group 
F43m). The end compound UCU5 possesses two crystallographically inequivalent 
copper sites in the ratio 4 : 1 at the 16e and 4c positions (Wyckoff notation), which 
are filled by Pd ions in the alloys. Therefore stoichiometric UCu4Pd could crystallize 
in an ordered structure, with only 4c sites occupied by Pd ions, in which case disorder 
might not play a role in the NFL behaviour of this compound. 

Chau et al [38] reported clastic neutron diffraction measurements on members of 
the \JG\x^- x Yd x series, and found no evidence for Pd/Cu disorder for x = 1. But the 
neutron scattering cross sections for Pd and Cu are similar, and the diffraction data 
do not rule out the possibility of interchange between Pd and Cu sites at the level of 
~4% occupation of 16e (Cu) sites by Pd atoms [and therefore ~16% occupation of 4c 
(Pd) sites by Cu atoms]. 

TF-/xSR linewidths have been measured in a number of samples of UCu4Pd [35], 
including a previously-studied powder sample [7] (Powder #1), the powder sample 
used for the neutron diffraction experiments (Powder #2), and a bulk polycrystal 
consisting of a few single crystals. Figure 2 plots 5K/(a*x), with a* calculated in the 
SRC limit, vs. the bulk susceptibility x, with temperature an implicit parameter [33]. 
(A combination of TF-/xSR and NMR results has shown that the SRC limit is 
appropriate to this system [7].) The data exhibit the near proportionality to x 
expected from the simple heuristic argument given above. (The dropoff for x ^ 
5 x 10~ 3 emu/mole (i.e., for T > 100 K) is not well understood but could be due 
to thermally-activated muon diffusion, which would motionally average the TF-/iSR 
line and reduce its width.) The samples were all prepared differently and are likely to 
have a range of defect concentrations, in which case the susceptibility inhomogeneity 
is remarkably insensitive to the amount of structural disorder. We note that the low- 
temperature values of 8K/(a*x) m UCu4Pd and UCu3.5Pdi.5 (which is definitely a 
disordered alloy) are also very similar [7]. Shown for comparison in figure 2 is Sx/x 
from the Kondo-disorder model [6-8] . The disorder-driven Griffiths-McCoy singularity 
theory (not shown) is also in good agreement with the data [12]. 

X-ray absorption fine structure (XAFS) studies in UCu4Pd have confirmed that 
Cu/Pd site interchange is substantial. Initially the observed disorder was argued to 
be sufficient to produce the required broad distribution of Tk in the single-ion Kondo 
disorder model [39], but subsequent annealing studies [40,41] call this into question. 
The XAFS results nevertheless leave little doubt that there is considerable interplay 
between site-exchange disorder and NFL behaviour in UCu4Pd. 
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Figure 2. Dependence of 6K/(a*x) on bulk susceptibility x, with temperature 
an implicit parameter, in UCu4Pd. Powder #1: previously-studied powder 
sample reference [7]. Powder #2: powder sample used in neutron diffraction 
studies reference [38]. Curve: &x/x from the Kondo-disorder theory reference 
[6-8]. From reference [33]. 



3. Muon spin relaxation and NFL spin dynamics 

3.1. Time-field scaling 

We have seen that disorder-driven mechanisms have been considered for the NFL 
properties of some f-electron systems [9,12], and it is natural to consider the 
possibility of extremely disordered or 'glassy' behaviour in the spin dynamics of these 
systems. On theoretical and experimental grounds it is believed that glassy dynamics 
lead to long-time correlations with distinct signatures as the freezing or 'glass' 
temperature T g is approached from above [19]. In a spin glass the spin autocorrelation 
function q(t) = (Si(t)-Si(0)) is theoretically predicted to exhibit power-law (q(t) — 
ct~ a ) or 'stretched-exponential' (q(t) = cexp[— (Xt)^]) behaviour [42]. Power-law 
correlation has been found in spin-glass AgMn using muon spin relaxation [19]. 

In a magnetic field H applied parallel to the muon spin direction muon spin 
relaxation (LF-/iSR) is due to thermally-excited f-electron spin fluctuations that 
couple to the muons. A muon at a given site experiences a time-varying local 
field -Bi oc (i) due to fluctuations of neighboring f moments. Following Keren et al 
[19,43], under motionally narrowed conditions the local muon asymmetry G(t,H) 
relaxes exponentially: 

G(t, H) = exp [-2A 2 r c (H)t] , (3.1) 

where A 2 = 7 2 ,(|i?i oc | 2 } is the time- averaged mean-square coupling constant in 
frequency units, 7^ is the muon gyromagnetic ratio, and the local correlation 
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time t c (H) is given by 



(H) = / dtq^co^u^t) = cu c (H) 



Jo 



(3.2) 



here = 7^-ff is the muon Zeeman frequency. We consider A and the prefactor c 
but not the functional form of u c (H) to vary from site to site in a disordered material. 
Then the sample-averaged asymmetry G(t, H) is given by 



where p(A, c) is the joint distribution function for A and c. A plot of In G(t, H) vs. t 
always exhibits upward curvature, because the initial slope on such a plot is the spatial 
average of the rate whereas the long-time slope is dominated by regions with low rates. 
Such upward curvature is difficult to achieve by other means, and is therefore evidence 
for spatial inhomogcncity in c and/or A. 

From equation (3.3) the time and field dependence of In G(t,H) enter in 
the combination u c (H)t. This means that G(t,H) scales as this combination 
independently of the form of p(A,c). For both the power-law and stretched- 
exponential forms of q(t) 



from equation (3.2) [in the high-field limit for stretched-exponential q(t)] [19]. Then 
u c (H)t scaling of G(t, H) results in the time-field scaling relation G(t, H) = G(t/H 7 ), 
i.e., a plot of G(t,H) versus t/H 1 is universal for the correct choice of 7. The sign 
of 7 — 1 distinguishes between power-law (7 < 1) and stretched-exponential (7 > 1) 
correlations [19]. It is important to note that no particular form for G(t) is assumed 
in this analysis. 

As noted above, the early- and late-time slopes of In G(t, H) are determined by 
spatially different sets of muons. If G(t) obeys time-field scaling, the shape of G(t) 
is independent of field, which implies that t c (H) [cf. equation (3.2)] has the same 
functional dependence for all sites in the sample. Thus observation of time-field 
scaling is evidence that the local dynamics exhibit scaling, with an exponent that 
is the same at all sites even if the sample is disordered. The significance of this 
property is discussed in section 5.1. 

3.2. UCu 5 - x Pd x , x = 1.0 and 1.5 

3.2.1. Muon spin relaxation [16]. Evidence from LF-/1SR experiments suggests that 
spin dynamics in the NFL alloys UCus-^Pda,, x = 1.0 and 1.5, are glassy, i.e., strongly 
slowed by disorder [44]. The sample average muon spin relaxation function G(t,H) 
(i.e., the asymmetry in muon counting rate; see, e.g., [45]) is strongly sub-exponential, 
indicating a quenched inhomogencous distribution of relaxation rates as discussed 
above, and obeys the time-field scaling relation G(t,H) — G(i/iJ 7 ) for applied 
magnetic field H between ~15 Oe and ~1 kOe. The field dependence corresponds 
to a measurement of the Fourier transform of q(t) over the frequency range w 
200 kHz-14 MHz, where 7^ = 2ir x 13.55 kHz/Oe is the muon gyromagnetic ratio. 
Power-law behaviour of q(t) is implied by the observation 7 < 1 [19], and also by the 
temperature-frequency scaling found in the inelastic neutron scattering (INS) cross 
section [46] (the LF-/1SR and INS results are compared in more detail in section 3.2.2). 
The LF-/iSR time-field scaling measurements extend the frequency range over which 




(3.3) 



Uc(JT) oc H-'< 



(3.4) 
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power-law correlations are observed in UCus^Pdz downward by three orders of 
magnitude. In addition, zero-field /zSR (ZF-juSR) shows no sign of static magnetism 
or spin freezing in UCi^Pd above 0.05 K. This together with glassy scaling points to 
some form of quantum spin-glass behaviour [21,47]. 

Figure 3 gives G(t) in UCu3.5Pdi.5 for T = 0.05 K and values of applied field H 



UCu 35 Pd 15 7= 0.05 K 




2 4 6 8 

Time (jxs) 

Figure 3. Field dependence of sample-averaged muon asymmetry relaxation 
function G(t) in UCu3.sPdi.5, T = 0.05 K. Curves: fits as described in text. 
From reference [16]. 

between and 2.5 kOe.f The relaxation slows with increasing field. For low enough 
fields we expect the field dependence to be due to the change of uo^ [cf. equations 
(3.1)-(3.3)] rather than a direct effect of field on q(t); a breakdown of scaling would 
occur for high fields where this might cease to be true. 

The same asymmetry data are plotted in figure 4 as a function of the scaling 
variable t/H 1 . For 7 = 0.7 ± 0.1 the data scale well over more than three orders of 
magnitude in t/H 1 and for all fields except 2.5 kOe. Fields hbH > k^T would be 
expected to affect the spin dynamics, and indeed the static susceptibility of UCu4Pd is 
suppressed by fields ~1 kOe below ^0.5 K (Vollmer et al [48]). The scaling exponent 7 
is less than 1, implying that q(t) is well approximated by a power law (or a cutoff 
power law [19]) rather than a stretched-exponential or exponential. From the /iSR 
data q(t) « ct~ 3±01 . We note again that no specific form for the muon asymmetry 
function has been assumed. A scaling plot is given in figure 5 for UCu4Pd, T = 0.05 K. 
Here the scaling exponent 7 = 0.35 ±0.1 is significantly smaller than in UCu3.sPdi.5. 
Scaling again breaks down for high enough fields. 

The fluctuation-dissipation theorem J relates t c (H) to the imaginary compo- 

f Figures 3—6 are reprinted with permission from D E MacLaughlin, O O Bernal, R H Hcffner, G 
J Nieuwenhuys, M S Rose, J E Sonier, B Andraka, R Chau, and M B Maple, Phys. Rev. Lett. 87, 
066502 (2001). Copyright © (2001) by the American Physical Society. 

I The fluctuation-dissipation theorem may be invalid in the frozen spin-glass state because of broken 
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Figure 4. Data of figure 3 plotted versus the scaling variable t/H 0,7 . From 
reference [16]. 
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Figure 5. Scaling plot of G(t, H) vs. t/H - 35 for UCu 4 Pd, T = 0.05 K. From 
reference [16]. 



crgodicity, but presumably holds in the high-tempcrature state considered here. See, e.g., Fischer K 
H and Hertz J A 1991 Spin Glasses (Cambridge: University Press) p 138. 
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nent x"( w ) 01 the local (g-independent) f-electron dynamic susceptibility: 




X "(u), T) as u)-tF(hLj/k B T) , 



(3.5) 



(3.6) 



with 7 = 0.33 and F(x) = tanh(x/1.2) for both UCu 4 Pd and UCu 3 . 5 Pdi. 5 . Using this 
form in the limit id, t c (H) obtained from equation (3.5) is independent of T and 
proportional to H~ 7 ; the latter is in accord with the LF-/^SR time- field scaling. The 
INS value of 7 agrees with LF-^SR data at T = 0.05 K for UCu 4 Pd (7 = 0.35), but 
not for UCu 3 . 5 Pdi.5 (7 = 0.7) (but see section 3.2.2). 

To go further one must fit the LF-^iSR data to an appropriate functional form 
for the asymmetry. On purely empirical grounds we use the stretched-exponcntial 



Values less than unity of the 'stretching power' K (not to be confused with the Knight 
shift) yield sub-exponential relaxation corresponding to a distribution of relaxation 
rates (see above and, e.g., reference [19]). This form is used because it characterizes 
an a priori unknown relaxation-rate distribution, and because A conforms with the 
general definition of the relaxation time 1/A as the time at which G(t) decays to 1/e of 
its initial value. For H — (ZF-^tSR) the data were fit to the product of equation (3.7) 
and the zero-field Kubo-Toyabe (K-T) function [50] characteristic of static relaxation 
by nuclear dipolar fields at muon sites. This form is expected when the muon 
local field has both static nuclear dipolar and dynamic f-moment contributions. A 
nuclear dipolar field —2.3 Oe was measured in both alloys for T 3> 1 K, where the 
contribution of U-moment fluctuations to the zero-field muon spin relaxation rate 
becomes negligible. Nonzero values of H were chosen large enough to decouple the 
muon spin relaxation from the nuclear dipolar field [50] leaving only the dynamic 
U-moment contribution, so relaxation data for these fields were fit to equation (3.7) 
without the K-T function. The curves in figure 3 are examples of these fits. 

Figure 6 gives A(T) and K(T) for UCu 4 Pd at three values of H. With decreasing 
temperature A increases slowly and saturates to a constant below 0.11-0.2 K. As 
noted previously the INS scaling predicts a temperature-independent relaxation rate, 
in mild disagreement with the observed weak temperature dependence of A(T). The 
stretching power K is 0.4-0.5 at 1 K, indicative of a broad distribution of relaxation 
rates [19], and increases slightly with decreasing temperature. Similar behaviour is 
exhibited by A(T) and K{T) in UCU3.5PCI1.5 (data not shown), with rates slower than 
in UCu 4 Pd by -30% at low fields and -100% at 100-300 Oe due to the larger scaling 
exponent 7. 

No anomaly is found in the ZF-^SR data at temperatures —0.1-0.2 K, where 
specific heat and ac susceptibility measurements (in different samples) suggest 
spin-glass-like freezing [48,51]. The muon-f-moment coupling in UCus-^Pd^ is 
predominantly dipolar [7] with a coupling field 0.55 ± 0.05 kOe//iB- Randomly- 
frozen moments of the order of 1 /xb/U ion would result in a ZF-/xSR relaxation 
rate —50 /is -1 , two orders of magnitude larger than the observed rate. This result 
places an upper bound of — 10~ 3 ^b/U ion on any frozen moment in UCu 4 Pd or 
UCua.sPdi.s. The discrepancy with the results of [48] might arise from differences 
in annealing conditions, or from a difference in sensitivity of the two techniques 
to superparamagnetic clusters or inclusions [44]. Such clusters could contribute 



G(t) = exp[-(At) K ] . 



(3.7) 
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significantly to the magnetic susceptibility, but might not occupy enough sample 
volume to be observed in ZF-^iSR experiments. 



3.2.2. Scaling exponent: temperature dependence [17]. Time-field scaled LF-/iSR 
asymmetry data for UCu 3 5 Pd! 5 and UCu4Pd are given in figures 7 and 8, 
respectively, for temperatures of 0.05 K and 0.5 K§. In UCu 3 . 5 Pdi.5 7 from muon 
spin relaxation exhibits a temperature dependence, becoming smaller at higher 
temperatures: 7(0.05 K) = 0.7 ± 0.1, 7(0.5 K) = 0.5 ± 0.1. This behaviour is 
reminiscent of spin-glass AgMn above the spin- freezing 'glass' temperature T g , where 
7 increases as T — > T g [19]. It suggests a monotonic temperature dependence of 7 
in UCu3.5Pdi.5 between the low-temperature muon spin relaxation results and the 
higher-temperature (>10 K) INS value 7 = 0.33 [49]. 

In UCu4Pd, on the other hand, 7 = 0.35 ±0.10 from muon spin relaxation at both 
0.05 K and 0.5 K (figure 8). This value is in agreement with INS experiments above 
10 K [49], which give 7 = 0.33 as in UCu 3 . 5 Pdi. 5 . The temperature independence of 

§ Figures 7—8 are reprinted from Physica B, vol 312-313, MacLaughlin D E, Hcffner R H, Bcrnal 
O O, Nieuwenhuys G J, Sonier J E, Rose M S, Chau R, Maple M B and Andraka B, "Slow spin 
dynamics in non-Fermi-liquid VC\i5^ x Pd x , x = 1.0 and 1.5," pp 453-455, Copyright © (2002), with 
permission from Elsevier. 
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7 in UCu4Pd suggests that the slow fluctuations are quantum rather than thermal in 
origin. In UCU3.5PCI1.5 the Pd concentration x is closer to the value x ~ 2 above which 
a spin-glass phase is observed [36], and the proximity of this phase may be reflected 
in the low-temperature increase of 7. 

3.3. CePtSh- x Ge x [18] 

CePtSii-zGez is a NFL system [52] with similarities to CeCu6_ x Au x , except that 
the residual resistivity of the latter alloy is somewhat lower. The end compound 
CePtSi is a heavy-fermion paramagnet. Ge doping expands the lattice and favors 
antiferromagnetism. A nonzero Neel temperature T/v appears for x > 0.1, which is 
therefore a candidate for a QCP. But magnetic susceptibility and 29 Si NMR studies 
of CePtSi and CePtSio.gGeo.i [15] show that both alloys exhibit strong disorder in the 
susceptibility above 2 K, which when extrapolated to lower temperatures accounts for 
much of the NFL specific heat. 

LF-/zSR experiments below 1 K in CePtSii-^Gea;, x — and 0.1, exhibit 
rapid relaxation and time-field scaling behaviour qualitatively similar to that seen 
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in UCus-ajPda;, x = 1.0 and 1.5 [18]. Figure 9 gives the muon decay asymmetry 
relaxation function for CePtSio.gGeo.i at T = 0.05 K and several field values. 
Normalized relaxation functions G(t, H) for CePtSii-^Gea; at T — 0.05 K are plotted 
in figure 10 against the scaling variable t/if 7 ||. For CePtSi the best scaling is achieved 
with 7 = 1.6 ± 0.1. As before, this procedure does not require knowledge of the 
functional form of G(t, H). For CePtSio.gGeo.i scaling is not obeyed for H < 300 Oe, 
but above this field scaling is found with the same value of 7. This value is considerably 
larger than in UCus^Pdz for T < 1 K, where 7(35=1.0) = 0.35 and 7(1=1.5) = 0.5- 
0.7 [16,17] (cf. section 3.2). 

The data are well described by the stretched exponential function Aexp[— (At) K ] 
introduced in section 3.2.1, as shown in figure 9 for CePtSio.gGeo.i- (For both Ge 
concentrations K is constant at about 0.25 for fields above ~50 Oe.) The field 
dependence of A is plotted in figure 11 for T — 0.05 K. For x — A varies as 
H r , 7 w 1.6^f, whereas for x — 0.1 this is true (with roughly the same exponent) 
only above a crossover field of about 100 Oe. This behaviour and the fact that 7 is 
greater than 1 are both consistent with a stretched-exponential spin autocorrelation 

I Figures 9-11 and 14 are reprinted from Physica B, vol 326, MacLaughlin D E, Rose M S, Young 
B L, Bernal O O, Heffner R H, Morris G D, Ishida K, Nieuwenhuys G J and Sonier J E, "/xSR and 
NMR in f-electron non-Fermi-liquid materials," pp 381-386, Copyright © (2003), with permission 
from Elsevier. 

f For G(t, H) = G[A(H)t] time-field scaling G(t, H) = G(t/H~t) implies A(H) oc ff-T. 
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Figure 9. Semi-log plots of LF-^tSR asymmetry data in CePtSio.gGeo.i, T = 
0.05 K. Curves: fits to stretched exponential Aexp[—(At) K ]. From reference [18]. 

function q(t) — exp[— (t/r c ) 7_1 ] [19, 20], rather than the power-law q(t) = t 1 ^ 1 (7 < 1) 
found in UCus-^Pda;. The crossover from A w const, to A a iJ~ 7 occurs for 
u^t c ~ 1. The data of figure 11 indicate that r c is considerably longer for CePtSi than 
for CePtSio.gGeo.i, due perhaps to less pinning of fluctuations in the more ordered 
end compound. This difference in behaviour is related to details of the microscopic 
mechanism for the disordered spin dynamics, and is not well understood. 

34. YbRh 2 Si 2 [53-55] 

The NFL compound YbRh 2 Si2 seems to be a suitable system for the study of 'clean' 
NFL physics. A number of bulk measurements [30] indicate that YbRli2Si2 exhibits 
antiferromagnctism with a very low Neel temperature Tn ~ 70 mK, with a field- 
induced QCP at the orientation-dependent critical field H c that suppresses Tn to 
zero. Resistivity and specific heat measurements at low temperatures show Ap oc T 
and C c \/T oc — InT over a temperature range of more than a decade. The NFL 
behaviour is suppressed and Fermi-liquid behaviour is recovered by the application of 
fields > H c . These results suggest that in zero field YbRh 2 Si2 is quite close to a QCP. 

Figure 12 gives the relaxation function G(t, H) for the muon decay asymmetry 
in YbRh 2 Si2 at T = 20 mK [53] + . At zero field the data can be fit to a static 
Kubo-Toyabe (K-T) relaxation function [50] corresponding to a random distribution 
of static local fields with rms width « 2 Oe. These static fields can be attributed 

+ Figures 12 and 13 are reprinted from Physica B, vol 326, Ishida K, MacLaughlin D E, Bcrnal O O, 
Hcffncr R H, Nicuwcnhuys G J, Trovarclli O, Geibel C and Steglich F 2003, "Spin dynamics in a 
structurally ordered non-Fermi liquid compound: YbRli2Si2," pp 403-405, Copyright © (2003), with 
permission from Elsevier. 
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to weak ordered Yb moments (10~ 3 -10~ 2 fis) in the antiferromagnetic state (Neel 
temperature Tn — 70 mK), since the K-T rate is much larger than expected from 
nuclear dipolar fields and sets in only below Tjv- The relaxation data for nonzero 
applied field can be fit to an exponential relaxation function G(t, H) = exp[— t/T\{H)]. 
An applied field of 11 Oe is more than five times larger than the estimated field at 
the muon site due to static Yb-moment magnetism [55], and therefore is large enough 
to decouple this static field. Thus the relaxation observed for H > 10 Oe is dynamic. 
Therefore the field dependence of the rate 1/Ti is not due to decoupling, but instead 
suggests a significant frequency dependence to the local-field fluctuation spectrum at 
the low muon frequencies. 

It is important to note that the observed exponential form for H > 11 Oe is 
evidence that the relaxation rate is substantially uniform throughout the sample, 
since as described above in section 3.1 the signature of relaxation-rate inhomogeneity 
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is a sub-exponential relaxation function, i.e., upward curvature of G(t) on a semi- 
log plot [16]. The muon spin relaxation is therefore probing spin fluctuations in a 
structurally ordered NFL system. 
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Figure 13 gives the dependence of 1/T\ on longitudinal field at T = 20 mK. The 
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zero-field rate given by the arrow is the static K-T value obtained from the fit shown 
in figure 12. The relaxation rate shows a weak field dependence for magnetic fields 
less than ~10 Oe but varies more strongly, as H^ 1 , 7 = 1.0 ± 0.1, for higher fields. 

As discussed in section 3.2.1, this behaviour is consistent with a scaling law of the 
form x"{lo,T) oc lo^ 1 F{ui/T), F(x) — > x for small x. This scaling leaves l/Ti(w,T) 
independent of temperature, in rough agreement with the observed weak (~ — InT) 
dependence at H = 19 Oe (not shown; see reference [54]). Furthermore, for 7 as 1 q(t) 
varies slowly with t [19]; the low-temperature spin fluctuations are very long lived. 

The exponential behaviour of the relaxation function in the present LF-/1SR 
experiments indicates that YbRh2Si2 is an ordered stoichiometric compound. Thus 
disorder-driven theories [9, 13] seem to be ruled out for this material. The /iSR results 
strongly suggest that YbRli2Si2 is a compound in which NFL behaviour is induced by 
homogeneous critical spin fluctuations. 

4. Effect of disorder on NFL behaviour in heavy-fermion metals 

We next consider relations between LF-/zSR relaxation and other properties in a 
number of NFL systems. Figure 14 gives the muon spin relaxation functions in 
CeNi2Ge2, YbRli2Si2, CePtSi, and UCu4Pd at low temperatures [18]. In these 
experiments a small longitudinal magnetic field has been applied, if necessary, to 
decouple any nuclear dipolar field, so that the relaxation is purely dynamic. The 
estimated relaxation functions for CeCu5.gAu .i [24] and Ce(Ruo.5Rh .5)2Si2 [26,27] 
are comparable to that for CeNi2Ge2, i.e., the relaxation rates are also very low in 
these systems. CePtSi and UCu4Pd show the most rapid and most noncxponcntial 
relaxation, whereas muons in CeCus.gAuo.i, Ce(Ruo.5Rho. 5)2812, CeNi2Ge2, and 
YbRh 2 Si2 all relax much more slowly. YbRh 2 Si2 exhibits weak exponential muon 
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spin relaxation at low temperatures and fields, with a rate 1/Ti that varies inversely 
with longitudinal field [53] (cf. section 3.4). 

Thus there seems to be a range of relaxation behaviour in NFL materials. Is this 
due to differences in characteristic energy scales, or is disorder intimately involved? 

The low-temperature muon spin relaxation data for these materials have been 
fit to the stretched-exponential form G(t,H) = exp[— (At) K ] as described above. 
Normalizing A by t^j, where u mo i is the volume per mole of f ions, roughly accounts 
for the concentration of f-ion relaxing centres*. Figure 15(a) gives plots of v^^A 
and K vs. the residual resistivity p(0) [18,32]. A smooth and marked increase of 
w^ ol A with increasing p(0) is observed, together with a decrease of K from ~1 to 
~l/2. There appears to be a crossover from ordered to disordered behaviour for 
p(fy ?t 100 cm ; disorder-driven effects are small even in solid solutions such as 
CeCus.gAuo.i [p(0) ~ 30 cm]. This suggests that disorder is important in, and 
could even be the driving mechanism for, slow f-electron NFL spin dynamics. 

We also wish to examine the question of correlation between the relaxation 
behaviour and some measure of the characteristic fluctuation rate, or characteristic 
energy scale, of the system. By definition NFL metals at low temperatures exhibit no 
characteristic energy other than the temperature. Nevertheless one can use a quantity 
such as the low-temperature Sommerfeld specific-heat coefficient 7(T) = C(T)/T (not 
to be confused with the time-field scaling exponent) as a rough gauge of the expected 
fluctuation rate, and hence of muon spin relaxation behaviour if slow fluctuations are 
correlated with large values of 7(T). In NFL materials 7(T) is temperature dependent, 

* The dipolar interaction between a muon and an f moment at distance r varies as r~ 3 oc u"^, as 
does the RKKY indirect interaction in simple models. This interaction enters the transition rate 
as the matrix element squared, so that the f-ion concentration dependence is normalized out after 
multiplication by f^ ol . Our conclusions arc not sensitive to this normalization. 
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and we have arbitrarily chosen to use the value 7(T=1K). The conclusions do not 
depend critically on this choice. 

In figure 15 (b) the relaxation parameters are plotted against j(T=l K), which is 
admittedly only a crude estimator of energy scales in these systems. Nevertheless, in 
a conventional picture smaller values of 7(T=1 K) would imply larger characteristic 
energies, faster spin fluctuations, more motional narrowing, and slower muon spin 
relaxation rates. This is not observed; there is no apparent correlation between muon 
spin relaxation parameters and -f(T=lK). It appears that disorder slows down the 
spin fluctuations, much as in the paramagnetic state of spin glasses but with no spin 
freezing down to the lowest temperatures of measurement (~20 mK). We appear to 
be dealing with spin glass behaviour in the extreme quantum limit [21,56]. 
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5. Concluding remarks 

5.1. Local or cooperative dynamics? 

The muon spin relaxation results in disordered systems can be compared with existing 
local (single-ion or cluster) theories of disorder-driven NFL behaviour. The Kondo- 
disorder model cannot account for the order of magnitude of the experimental muon 
spin relaxation rates at low temperatures [33]; there simply are not enough slowly- 
fluctuating uncompensated moments at low temperatures to account for the observed 
muon relaxation rates. The situation for the Griffiths-McCoy singularity model is 
more complicated, however. 

In an early version of the theory [12], which treats the effect of f- moment 
clustering, there is no dissipation in the f-electron spin dynamics, and the local cluster 
dynamic susceptibility is sharply resonant at a distributed characteristic tunneling 
energy E: 



Note that this form does not exhibit frequency scaling x"( w ) <x uo^ 1 F(uo /T). Together 
with the distribution function P(E) oc _E~ 1+A , where A < 1 is a nonuniversal scaling 
exponent, we have 



This scaling form is in agreement with the INS data for A w 0.7 and is reminiscent of 
the LF-^iSR results. 

But in LF-/iSR experiments it is the muon spin relaxation function G(t, H), not 
x"{oj), that is spatially averaged. As discussed in section 3.1, the observed time-field 
scaling of the spatially-averaged muon spin relaxation function G(t, H) shows that 
the local x"{un), not just the spatial average x"( w ); scales as w~ 7 . Local frequency 
scaling is not a property of the Griffiths- McCoy singularity model, in which the scaling 
is found only after the average has been taken, and only for a sharply resonant (and 
nonscaling) form [equation (5.1)] of the local x"( w : E). A later form of this theory [13] 
considers dissipative effects, which broaden x"( w > E) but do not give it a scaling form. 
Furthermore, if the width of x"( w ) is much greater than u it is not hard to show 
that x"(uj) no longer follows P(u}—E), so that this mechanism for scaling of x"( w ) is 
lost. It is also difficult to see how the dynamic susceptibility of [13] would yield the 
observed temperature dependence of A(w Jli ,T). Thus the Griffiths-McCoy singularity 
theory does not seem to account for the /xSR results for a number of reasons. 

More generally, disorder-driven pictures that ascribes NFL behaviour to 'rare' 
objects (low-Tk spins or large clusters) lead to spatial distributions of spin fluctuation 
rates (slow for the rare objects), and thus are not in agreement with time-field 
scaling. We conclude that local disorder-driven models account for inhomogeneous 
distributions of static susceptibility in disordered NFL materials (section 2), but that 
the distributed local dynamics in these models do not describe the observed muon 
spin relaxation. Time-field scaling of the relaxation data is strong evidence for a 
cooperative mechanism, i.e., spin dynamics with a single correlation function for the 
entire sample. 



X "(u), E) cx S(uj - E) tanh(£/2T) . 



(5.1) 




(5.2) 
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5.2. Glassy behaviour? 

The question of local vs. cooperative dynamics seems to be answered in favor of the 
latter as discussed above, at least for the UCu5_ x Pd x system which has been studied 
most intensively. In contrast, a number of more general (and related) questions remain 
unanswered, or at best only partially answered: 

• Disorder slows down the (cooperative) spin fluctuations in disordered NFL 
systems [figure 15(a)]. Are these slow 'glassy' fluctuations ordered-system 
quantum critical fluctuations strongly modified by disorder, or new 'quantum 
glass' excitations created by disorder? Is there even a distinction in principle 
between these notions? 

• The correlation shown in figure 15 (a) strongly suggests that glassy dynamics 
depend 'universally' on residual resistivity. What does this mean physically? Is 
there a picture in which disorder-driven slowing down of fluctuations in a metal 
depends only on the electronic mean free path? In this regard the increase of 
impurity scattering near a QCP predicted by Varma, Miyake, and co-workers [57- 
59] may be relevant if there is a slowing effect of this scattering on spin lifetime 
in the system near criticality. 

• More generally, how are the fluctuations slowed down? Is frustration involved, 
in the sense that it is involved in slow spin fluctuations in spin glasses and 
geometrically frustrated systems? 

• Is there a mixing of (quantum) critical and Griffiths-McCoy singular behaviour? 
Do we need to invoke aspects of both phenomena [60]? 

The static and dynamic response of spins in f-electron NFL materials exhibits 
very strong effects of disorder, revealing poorly-understood instabilities of quantum 
criticality in these systems to strong inhomogeneity and glassy behaviour. It is to be 
hoped that the present experimental results lead to further endeavours in this area. 
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